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Abstract. The DA white dwarfs in the Sloan Digital Sky Survey, as analyzed in the papers for 
Data Releases 1 and 4, show an increase in surface gravity towards lower effective temperatures 
below 11500 K. We study the various possible explanations of this effect, from a real increase 
of the masses to uncertainties or deficiencies of the atmospheric models. No definite answer is 
found but the tentative conclusion is that it is most likely the current description of convection 
in the framework of the mixing- length approximation, which leads to this effect. 



1. Introduction 

Almost every contribution at this white dwarf meeting mentions the Sloan Digital Sky Survey 
(SDSS) or is even based on its data. Because of the huge impact the survey has on this field, it 
is very important to critically analyze all aspects of the data collection and reduction process, to 
understand possible uncertainties and biases. One problem rel ated to white dwarfs, w hich was 



obvious already in the first white dwarf catalogue publication ( Kleinman et al. . 20041 ) . was the 



apparent systematic increase of the mean surface gravity of DAs below effective temperatures of 
approximately 12000 K, determined from spect ral fits, was again demonstrated i n the data 
of the Data Release 4 ( Eisenstein et al. . 20061 ) and confirmed by Kepler et al. ( 200?! ) and 



DeGennaro et al.l (|2008l ). 



While the huge number of objects and the homogeneity of the survey make this feature 
very obvious, it is by no means a ne w discovery. A similar trend was firs t recognized with 
data from the Palomar-Green Survey ( Bergeron et al. ^ l990l: BergeronUl992l. and several other 
papers of the Montreal group), and more recently by lLiebert e t al. (2005). A si milar effect was 
also found in the data from the ESO Supernova la Progenitor Survey (SPY) by Voss ( 20061 ). A 
first conclusion therefore is that this effect is not a peculiarity of the SDSS data collection or 
reduction processes. In this paper we study the remaining possibilities 

• the effect could be real, i.e. the DAs show a real increase of the average mass towards lower 
temperatures 

• it could be caused by the fitting between observations and theoretical models used to 
determine the atmospheric parameters 

• uncertainties in input physics or a flaw in the models used 



For our study we use the data of D R4 as provided from th e SDSS database, with shght changes 
in reduction software compared to Eisenstein et al. ( 20061 ). We also use very nearly the same 



models as used in that work for our "standard analysis" (see below). The major difference are 
the fitting algorithms: whereas the SDSS uses the whole spectra, adjusted with a polynomial 
using information from the photometry, we here use only the spectral lines, with the continuum 
adjusted on both sides of the line center. This technique minimizes the influence of uncertainties 
in the flux calibration. The SDSS photometry gives an independent estimate of parameters, 
which is used to distinguish between the two solutions, which exist for most objects between 
8000 and 18000 K. 

In the main part of this study we vary some of the parameters in the input physics, in the 
hope of finding some hints for the origin of the effect. This concerns the treatment of convection, 
of the non-ideal effects in the occupation numbers, and a possible admixture of invisible helium. 
All these possible explanations have been proposed and studied in various papers of the Montreal 
group. We partially repeat this work here with a much larger and very homogeneous sample. 



2. Statement of the problem 

Figure [T] shows a subsample of the brighter, cool (< 16000 K) DAs ir i the DR4 result s . The 
parameters were determined by SJK; the only significant difference to Eisenstein et al. ( 20061 ) 
was a larger model grid. While the scatter of the distribution and trend towards large logs' is 
similar for the two subsets, it is much smaller than in the total sample, which includes many 
fainter stars (see DR4 for a comparison). A S/N of at least 20, and preferably > 30 is a minimum 
requirement for reliable spectral analysis; we will therefore in the remaining tests use the bright 
sample with g < 19. The trend in log g becomes even clearer, if we bin the results in Tefj intervals 
of 500 to 1000 K (Fig. [2]). The spectroscopic fits show a clear increase below Tefj = 11500 K, 
which is somewhat more pronounced in the DR4 results than in our own fits. Since everything 
else is very similar, we attribute this difference to the different fitting methods. A tentative 
explanation could be that the SDSS polynomial correction to the flux calibration is problematic 
for the very extended wings of the broad Balmer lines, especially in this temperature range. 
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Figure 1. Distribution of logg with T^s for 
a subsample of cool DAs from DR4. Crosses 
(red): 347 DAs with magnitude g < 18; circles 
(black): 360 DAs with g between 18 and 19. 



Figure 2. Average logs' values and 
distribution in bins of Tgff. Squares (black): 
original results from DR4; circles (red): our 
new fits to the spectra; horizontal bars (blue): 
results from photometry for the same objects. 
The vertical lines in all cases are the intervals 
containing the "inner" 50% of the distribution, 
cutting 25% on either side. 



The surface gravity derived from fitting the 5 SDSS magnitudes simultaneously are slightly 
lower than the spectroscopic results. They also show a small increase of log 5 for lower Teg 
but this is much less pronounced than in the spectroscopic data. We would like to make a few 
remarks regarding the SDSS photometry, which is not used as much for the analysis of white 
dwarfs as it deserves. In the case of the DAs the discrimination for T^s and log 5 is very good in 
the range from 16000 K to 8000 K, and a determination using e.g. the two-color diagram u — g 
vs. g — r can in principle give results as accurate as from spectroscopy, without the problem 
of having two possible solutions above and below the Balmer line maximum. However, there is 
still some small uncertainty in the calibration of the magnitude zero-points. As a test, we have 
increased the observed g magnitude by 0.03, which results in a perfectly constant average log 5 
sequence close to 8.0. While we do not imply that the g uncertainty is as large as this, we believe 
that because of the calibration uncertainties the photometry should at this time not have very 
high weight in a study of systematic effects in the stellar parameters. We will, however, in the 
following sections also show the effect of model changes on the photometric solutions, which in 
an ideal case would agree with the spectroscopic solution. 



3. Could the mass change v^iih Teg be real? 

The input physics for our "standard" models is the same as used for the DR4 fits and will be 
described in more detail below. Our spectroscopic fits with these models to the ~ 700 DAs 
brighter than g = 19 and Tefj between 8500 and 16000 K is our reference set for this study. 
There are 61 DAs with multiple spectra in this set; from the differences in the determined 
parameters we estimate a typical uncertainty of 2.4% for Teg, and 0.15 in \ogg. Since the 
instrumentation, reduction, models, and analysis methods were the same, this is certainly still 
an underestimate of the real errors. Nevertheless, the number of objects are large, and the log 5 
increase so systematic that it certainly cannot be explained by uncertainties of the parameter 
determinations . 

With our reference fits the mean mass of the DAs is 0.573 Mq for the objects above 11500 K 
and 0.658 below. Could this increase of the white dwarf masses be a real effect? Possible 
explanations might be different cooling times for high versus low mass white dwarfs, or higher 
mass progenitors (and consequently high mass wds) for older stars, which originated from older 
progenitor s with lower inetal ab undances. The latter argument is, however, contradicted by the 
findings of [Hansen et al.l (|2007l ). that old white dwarf s in globular c luster s have masses close to 
0.5 Mq. All these arguments have been discussed in iKeoler et all (|2007l ^ and found not to be 
convincing. It is also very difficult to explain the rather sudden increase around 11500 K with 
these explanations, which all should lead to a more gradual transition if any. 



4. Non-ideal effects in equation of state and occupation numbers 

The standard treatment of these effects in white d warf atmosphere models i s the Hummer- 
Miha l as-Dappen occupatio n probability mechanism (jHummer &: Mihalad . Il988l : iMihalas et al ' 



19881 : Dappen et al. . 1988). For any ef fect above 10000 K only the interactions with charged 



particles are of interest. Bergeroiil ( 19931 ) has demonstrated that a more consistent fit between low 
and higher Balmer lines can be achieved by increasing the parameter j3crit in the HMD formalism 
from 1 to 2. The effect is to decrease the absorption in the wings of the Balmer lines, compared 
to the standard treatment. We have experimented with this parameter and show two extreme 
cases in Figs. O and [H For the smaller value (stronger line wings) the spectroscopic solution 
results in smaller log 5, as expected, and the values agree better with the photometric solution. 
The opposite is true for larger f3crit- Neither of the changes offers any real improvements, since 
the jump at 11500 K remains. This could also be expected, because the effect of changing the 
strength of the charged particle interaction should be continuous around this critical region. We 
therefore exclude this change as a possible solution. 
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Figures. Average of log g as a function of Teff Figure 4. Similar to Fig. [31 but PcrU 
using models with Pcrit decreased to 1. Circles increased to 5. 
(red): spectroscopic solution; horizontal bars 
(blue): photometric solution 



5. Invisible helium in the atmosphere 

This possible explanation was also first proposed by the Montreal group (iBergeron et alJ . ll99nl ). 



Why could this solve the problem of increasing log 5? The surface gravity, which along with 
Teff are the stellar parameters determined in an analysis, is only a proxy. The emission from a 
hot plasma is - in addition to the elemental abundances - determined by the thermodynamic 
variables, e.g. pressure and temperature. These variables are of course not constant throughout 
the visible atmosphere, but they can be described by characteristic values. For the temperature 
this is the effective temperature. The pressure is determined by the hydrostatic equation ^ = ^ 
Assuming a constant absorption coefficient (per mass) k this can be integrated to r = 1 to give 
P ~ g/K, showing that the pressure is inversely proportional to the absorption coefficient per 
mass. The latter is much smaller for helium than for hydrogen around 10000 - 12000 K, and 
therefore a significant helium pollution leads to higher pressure, mimicking a higher surface 
gravity. 

A possible source of helium could be helium dredge-up by the convection zone (CVZ), which 
(with the standard description of convection in white dwarfs) increases significantly in depth 
for effective temperatures below 11500 K. If the total amount of hydrogen in the outer layer 
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Figure 5. Size of the convection zone in DA 
white dwarfs function of Teg. The two 
continuous (red) lines are the upper and lower 
boundary, the gray area is the convection zone. 
The vertical axis is a mass scale, giving the 
logarithm of the fractional mass counted from 
the surface of the star. Note the sharp increase 
in depth around 11500 K. The mixing- length 
parameters for this calculation are ML2/0.6 
(see text) 



is less than about 10 of the total mass, such a mixing could happen near the interesting 
temperature range. 
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Figure 6. Effect of adding helium to the Figure 7. Adding hehum with an abundance 
atmosphere with a number ratio of He:H = of He:H = 1 by numbers. 
0.1. Symbols and colors are the same as in the 
previous figures. 



For simplicity we show in Figs. [6] and [7] the effect of adding a helium contamination with He:H 
= 0.1 and 1 by numbers in all models with Teff < 11500 K. Obviously this goes in the direction 
of a constant mass for all temperatures. In reality we would expect a helium pollution, which 
increases with decreasing Tefj, with an even more convincing result (removing the dip between 
11000 and 12000 K). 

Could such a helium contamination be detected directly? We remind the reader of 
two recently discovered objects, which provide some guidance - HS 0146+1847 and GD362. 
Both stars were at f irst classified as massive DAZ, with very stro ng Balmer lines and Call 
( Koester et al. . 2005 : Gianninas et al. . 20041 : Zuckerman et al. . 200?! ). Yet, at closer inspection, 
they turned out to be in fact helium-rich, with normal masses. In both cases high resolution and 
high S/N spectra were necessary to discover the weak helium lines, although the helium content 
is far larger than necessary to explain the average log 5 increase. In order to find more similar 
objects we have systematically searched ~ 1600 spectra from the SPY survey for weak helium 
lines, with no new discoveries. However, an abundance ratio of He:H = 1 or smaller would in 
general not be detectable at the typical S/N of that survey. We also note, that both helium-rich 
objects also have metal traces, indicating that some kind of accretion has happened and the 
stars are peculiar. Nevertheless, from the observational point alone, helium pollution could be 
the solution. 

From the theoretical side this w ould be very difficul t to understand. The latest 
asteroseismological study of ZZ Cetis ( Castanheira fc Keplei . 20081 ) finds a few objects with 
a hydrogen layer log M/^/M < —8 (none < —10), but the vast majority have much thicker 
layers where no mixing is expected in the ZZ Ceti range. 



6. The treatment of convection 

The description of the energy flux by convection is one of the few remaining problems in stellar 
astrophysics. In the field of white dwarfs it has become customary to increase the freedom in 
the mixing-length approximation with 3 more parameters in addition to the mixing-length itself. 
The different versions are distinguished as MLl, ML2, ML 3 , with the mixing length added as 
in MLl/a = 2.0 or short MLl/2.0 (see e.g. iFontaine et"all . Il98ll : iKoester et al.l . Il994i . for the 
nomenclature). For our current study the main point is that these versions differ in the efficiency 
of energy transport: the more efficient a version is the smaller can the temperature gradient be 
which transports the necessary flux. This temperature gradient in turn has a strong influence 



on the line profiles, and to a lesser degree also on the continuum energy distribution. We have 
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Figure 8. Convection completely suppressed, 
energy transport by radiation only. 



Figure 9. Temperature gradient forced to 
adiabatic in convective regions. 



experimented with various mixing-lengths in the framework of the MLl and ML2 versions. The 
changes always bring changes to the sequence of average log (7, but none solves the problem in 
the sense that a constant average results. As two extreme examples (which are outside the ML 
framework) we present in Figs. [8] and [9] a model set with convection totally suppressed (extremely 
inefficient) and another one with the convective gradient forced to the adiabatic gradient (infinite 
efficiency). In the second case this has an influence even at rather high temperatures, because 
there are thin convection zones, which are, however, normally inefficient and don't change the 
radiative gradient. 

While this result is not encouraging, we nevertheless want to argue, that incorrect description 
of the convective gradient is still the most likely and in our opinion the only possibility to account 
for the observed effect. 

First, there is a very striking coincidence near ~ 11500 K between the steep increase in the 
thickness of the convection zone and the increase in the average logs'. Of the other possible 
solutions only the helium contamination would be able to explain this coincidence. 

Second, and even more important: we definitely know that our current description of 
convection is incorrect. The current choice of parameters in the n iixing-length a p proxim ation 
(ML2/0.6, or very similarly MLl/1.75) was introduced by Koester et al.l ( 1994 ) and 
Bergeron et al. (|l995l l to allow a consistent parameter determination from optical and ultraviolet 
spectra. The effective temperatures found with this version for ZZ Cetis are in many cases so hot, 
that the thermal time scale of the CVZ is at most a few seconds, much shorter than the pulsation 
periods. On theo retical grounds, these timescales are expected to be of similar magnitude, 
and consequently iBerg eron et al.l (|l99,^ ) concluded that no single choice of ML parameters 
can simultaneously reproduce optical and UV spectra and the location of the instability strip. 
Convection must be more efficient in deeper layers than in the atmosphere, at least when studied 
within the ML framework. 

This same conclusion was reached completely independently by iLudwig et al. When 
comparing the temperature profile from a two-dimensional hydrodynamic simulation of a white 
dwarf with Tefj = 12600 K, log (7 = 8 with mixing-length models of the atmosphere and envelope 
they found that the mixing-length must increase with depth. 



7. Conclusions 

The frustrating result of this study is that we cannot propose any change of approximations 
in the input physics, which would clearly lead to the desired result - a constant average mass 



for all white dwarfs for all effective temperatures. Variations of parameters in the Hummer- 
Mihalas-Dappen occupation probabilities can be excluded. Helium pollution cannot be definitely 
excluded, but could be tested with high resolution, high S/N spectra of (apparent) DAs in the 
range 10000 - 12000 K. Because of the asteroseismology results discussed above we consider this 
solution unlikely. 

That leaves by default convection as the only suspect. Given the current amount of 
uncertainty about the correct description, and the contradictions in the current applications, we 
consider it quite possible that a more physically sound description of the convective flux and 
gradient might solve the problem. 

At the end of this study a final word of caution: fitting observed spectra with theoretical 
models is as much an art as an objective science. One should be very sceptical of results from 
spectra with S/N < 30, as well as DA spectra that do not include the higher Balmer lines with 
high S/N. Never believe the errors so readily calculated by a routine; these are errors which 
consider only statistical errors in the observed data. With modern CCD spectra from large 
telescopes systematic effects of reduction, the fitting process, and model deficiencies result in 
much larger errors. These can only be estimated in a meaningful way from comparing results 
for the same objects from different observations, and if possible, from different authors, with 
different fitting routines and model grids. 
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